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After a derivation of low-energy limit of QCD, being this a non-local Nambu-Jona-Lasinio model, we are able to 
show that confinement emerges as a two-loop correction to the gluon propagator. One-gluon exchange is not enough 
as recently shown in literature about studies on the gluon propagator in the Landau gauge. 
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1. Introduction 

Quantum field theory can be approached with ana- 
lytical techniques if a set of classical solutions to field 
theories are known. This can be understood when we 
look at standard techniques like small perturbation the- 
ory where one generally starts from classical solutions 
of the free theory. Here we show that such a set of solu- 
tions exists for a class of nonlinear field theories giving 
us the chance to analyze their behavior in the infrared 
limit. One of these is Yang-Mills theory. 

When the low-energy behavior of Yang-Mills theory 
is known one can approach quantum chromodynamics 
(QCD) in the same limit. Things stay this way as a 
knowledge of the propagator of the gluon can yield the 
proper low-energy limit of QCDflJSl- In this way, a 
non-local Nambu-Jona-Lasinio model is recovered that 
has already discussed in [3]. But, in our case, the 
form factor is completely fixed by QCD rather than phe- 
nomenologically guessed. Such a form factor agrees in 
a strikingly good way with the one obtained for an in- 
stanton liquid [4J]. This gives a sound evidence that in- 
stantons could play a relevant role for the ground state 
of the theory. Finally, it is essential to point out that 
another approach provides a similar link between QCD 
and Nambu-Jona-Lasinio model Q]. 

We will show in this paper how confinement can 
emerge in such a scenario, confirming a recent finding 
in literature that one-gluon exchange is not enough to 
grant it ffj, given a proper form for the gluon propaga- 
tor in the infrared limit and higher loops are considered. 



2. Classical fields 

The simplest theory to analyze, starting from a clas- 
sical standpoint, is the following 



op + = j. 



(i) 
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When j = an exact solution exists and is <p — 
fi (2 1 A)* sn(p ■ x + 0, i) being sn an elliptic Jacobi func- 
tion and n and 9 two integration constants. This so- 
lution holds provided the dispersion relation p 1 = 
/u 2 ^[A]2 holds. This appears to be a free massive so- 
lution notwithstanding we started from a massless the- 
ory. Mass originates from the nonlinear term but the 
coupling must be taken finite as a small perturbation ex- 
pansion fails to recover it. We want to solve eq. (Q} 
when A — > oo. In order to achieve this aim we consider 
a technique firstly devised in the '80s [7]. One considers 
(f> a functional of the current j and makes a power series 
of it. It is not difficult to realize that this argument is 
consistent when 

+ //A j d 4 x'd 4 x"A(x - x')[A(x - x")fj(x) + 0(j(x) 3 ) 

being A(x - x') a solution to the nonlinear equation 
□A(x - x') + A[A(x - x')] 3 = ^ x 6\x - x'). Such a 
series has an interesting leading term. This term, when 
introduced in a quantum field theory, will yield a Gaus- 
siam term in the generating functional that is the hall- 
mark of a trivial theory. We realize that this set of clas- 
sical solutions are peculiar to a trivial theory in a lead- 
ing order. This means that such a theory is manage- 
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able notwithstanding the presence of strong nonlineari- 
ties. But this program can only be complete if we are 
able to get the Green function of the theory. In this 
case this can be easily obtained by starting from the 
theory at d = 1+0. So, the equation to be solved is 
^A (f-f') + /l[Ao(f-f')] 3 = P- 2 S(t-t') and this is easily 
obtained (8J]. A Lorentz boost or a gradient expansion 
provides the complete solution 



A(p) = y ( 2„ + i)4-tl ) 



1 



n=0 



K 2 (i) 1 + e-(2n+l> p 2- m 2 + i e 



(2) 

being m„ = (2n + V)(n/2K(i)) (A/2)* (1 and K(i) * 
1.3111028777 an elliptic integral. This Green function 
appears to be consistent with a strong coupling expan- 
sion. 

The approach of current expansion devised in the ' 80s 
10] was conceived to manage QCD. The leading order 
for the Yang-Mills part of the theory was Gaussian. So, 
we try to see if, assuming to exist an infrared trivial fixed 
point for Yang-Mills theory, a set of classical solutions 
for Yang-Mills equations exists to support this approach 
as already happens for the scalar field. We start from the 
equations of motion 

d^d^A" - (l - ij dA&Ap + gf abc A h »(d,A c v - S v Ap 



+ gf abc ^(A b A c v ) 
+ if ahc f de A btl A d u A e v 



-Jv 



In this case, we note that the homogeneous equations 
can be solved by looking at instanton like solutions 
A^(x) = rf^<p(x) being rf'^ a set of constants. The problem 
here is that we have to fix the gauge. We consider the 
case of the Lorenz gauge that is equivalent to the Landau 
gauge in quantum field theory. In this case the equation 
collapses to d^d^ip + Ng 2 //) 3 = and our problem is 
equivalent to the one of the scalar field we discussed be- 
fore (we just note that this is not true for other gauge 
choices where we can only claim an asymptotic corre- 
spondence |0]). Then, we are able to write down im- 
mediately the gluon propagator in the Landau gauge, 
given eq.©, just setting A - Ng 2 and introducing the 
factor 6 a b (%>> ~~ PiiPvl p 2 )- This conclusion implies that 
instantons are relevant in the ground state of the theory. 

Our fundamental starting hypothesis, that Yang-Mills 
theory admits an infrared trivial fixed point, should be 
tested on the lattice to get support. This is indeed the 
case. From the German group ifioll the running cou- 
pling is seen to go to zero as momenta lower (see fig. 
[TJ. Similarly, the French group obtained such a behav- 
ior for the running coupling with a different definition 
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Figure 1: Running coupling for 64 4 and 80 4 atff = 5.7 full. 



Ill 10 . They show a perfect consistency with an instan- 
ton liquid model in agreement with the scenario we are 
depicting here. 

3. Quantum fields 

Quantum theory for the scalar field is man- 
aged with the generating functional Z[j] = 
N J [#] exp [i J d 4 x [(dcf>) 2 - Acf> 4 /4-]]. We rescale 

the space-time coordinates as x — > y/Ax and take a 
strong coupling expansion (f> = J^L A ~"(f>„. Then, 
we see that the leading order reduces to solve the 
equation D0o + A<p^ = j that we already know how 
to manage. This means that the leading order is just 
a Gaussian functional with the propagator given by 
eq.(f2| and the current expansion discussed above lfl2ll . 
Next-to-leading order can be also computed. Our 
fundamental result is that the massless scalar field 



theory in four dimensions is infrared trivial [12]. Mass 
spectrum is given by m„ = (2n + l)(n/2K(i)) (A/2)* fj. 
that represents free particles with a superimposed 
spectrum of a harmonic oscillator. For the Yang-Mills 
generating functional a similar Gaussian functional is 
obtained 



Zo[j] = N exp 



d'x'd 4 



'r(x')D ab ( X ' 



')/V 



) 



given the current expansion A^ — A J d A x'D ab Y (x - 
x')j bv (x') + 0(l/^Ng) + O(f) and the propagator 
D ab (p) = 6 ab (r, flv -^)A(p) being A(p) given by 
eq.(f2| with A = Ng 2 . The spectrum is given by free 
glue particles with a superimposed spectrum of a har- 
monic oscillator. Similarly, for the ghost field, our in- 
stanton solutions decouple it from the gauge field yield- 
ing a free particle propagator for a massless field. These 
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properties of the quantum Yang-Mills field represent the 
so-called "decoupling solution" II 1 31 - 4 1 511 . This solution 
is the one recovered in lattice computations |U 18]. 
Numerical solution of Dyson-Schwinger equations as 
given in IU3I1 (see fig© are in perfect agreement with 
it. 



v__ 



compare this form factor both with the oneguessed 
in y|] and the one from an instanton liquid [4] in fig. 
[3] The result is strikingly good for the latter show- 
ing how consistently our technique represents Yang- 
Mills theory through instantons. Finally, the cr field in- 




Figure 2: Comparison of our propagators with numerical solution of Dyson- 
Schwinger equations as in jj3:]. Running coupling is given on the right. 



Once we know the behavior of Yang-Mills theory, we 
are able to get the low-energy behavior of QCD. One 
can show that it takes the form [ 1 ] 



S 



= J d 4 x 



ta \2 1 2 2 

■jiocf) - -m cr 



(3) 



where the cr field arises from the gluon propagator in the 
Gaussian generating functional of the Yang-Mills action 
and is the contribution from the mass gap of the theory, 
being mo = (n/2K(i)) and & is the string tension 
O (440 MeV) 2 ). For the quark fields one gets 

S q = ^ I d 4 xq(x) \i$ - m q ] q(x) 
q J 



■f d 4 x' 



A(x - x')x 



2j 2j ^Wt^^MtT/'?'^)?^) 



A" 



(4) 



This result recovers the non-local Nambu-Jona-Lasinio 
model given in [ 3] but directly from QCD provided the 
form factor is 



0{p) 



1 



-g 2 A(p) = -C(p) 



(5) 



being A(p) that in eq.©, C(0) = 1 and 2£(0) = G 
the well-known Nambu-Jona-Lasinio coupling. So, the 
value of G is fixed by the gluon propagator. We can 



Figure 3: Comparison of our form factor with those provided in Hell et al.| 3] 
and instanton liquid 4] as given in [1]. 

traduced by the gluon field contributes to the cr field 
arising after bosonization in the Nambu-Jona-Lasinio 
model lfl9ll giving a renormalized coupling constant 
G e ff = 1/(1 /G + OTq). This field is effectively inter- 
acting with the other fields when these are given as 
bound states. This field can be identified with the cr 
or f0(600) resonance, whose mass has been estimated 
to nia- = 441+' 6 MeV S and 457+|^ MeV [21]. Our 
estimation, in the Nambu-Jona-Lasinio contact approxi- 
mation that provides m* =214 MeV for quark effective 
mass and m„ = 139.7 MeV, is = sj4m* 2 + m 2 = 
451 + 20 MeV. The error should be ascribed to the 
string tension that is estimated with an error of 20 MeV. 
Notwithstanding the rough approximation, this model 
provides a very good value for this mass. 

4. Confinement 

In order to verify if, at leading order (one-gluon ex- 
change), the scenario we depicted yields a confining 
theory we need to evaluate the integral for the interquark 
potential |6|] 



2C C 

V(x) = J d 3 pa s (p)Mp)e ipx . 



(6) 



Here, one has to use the running coupling ff s (p) that, in 
the infrared limit, we have proved to go like p 4 II lLllZl . 
This implies for the potential 



V(r) 



g.(0) 1 d 4 
A 4 7cM 



{-\) n e 



K 2 (i) 1 + e -(2»+i> 



(7) 
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and, due to massive excitations, one gets a screened po- 
tential. This appears to agree very well with the con- 
clusions given in l6[] but not in agreement with Cornell 
potential observed on the lattice for quenched simula- 
tions. So, if one just stops at the decoupling solution 
in the deep infrared, confinement is not recovered. This 
can be obtained if one is able to compute the next-to- 
leading order correction to the decoupling solution. For 
our scenario, this is presented in [22] yielding a two- 
loop sunrise diagram. The correction to the leading or- 
der propagator takes the form 



2 ri f d A p\ d 4 p2 B m 
F(A,p) = 2j J (27r) 4 (2w) 4 2 2 



2 p\-m 2 n2 



B ni 



(p-pi - pi) 1 -ml 



fixed point. A low-energy limit of QCD is so obtained 
that reduces to a non-local Nambu-Jona-Lasinio model 
with all the parameters and the form factor properly 
fixed by QCD. <x meson gets a value of mass in close 
agreement with recent studies on 7T7r-scattering and its 
very nature is mostly gluonic. Confinement can emerge 
in all this picture at two-loop level. Overall agreement 
with numerical studies of truncated Dyson-Schwinger 
equations and lattice studies make this scenario really 
compelling. 
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with m n . = {In + l)(n/2K(i)) 01/2)' A, A = Ng 2 and B n . 
can be obtained from eq.©. Techniques to evaluate this 
integral at small momenta, as we need, are well-known 
112 311 and the result can be stated into the form 



A 2L (p) = AO) 



(8) 

being c,- some inessential numerical constants. This 
correction provides the needed p 4 Gribov contribution 
to the propagator to get a linear term in the potential, 
taking into account the behavior of the running cou- 
pling that increases as the fourth power of momenta. 
Our scenario is in striking agreement with the numer- 
ical solution of truncated Dyson-Schwinger equations 
as shown in fig. [2] Truncated Dyson-Schwinger equa- 
tions where considered in [6] and our conclusions agree 
at one-gluon exchange level as they should. 

We note that a regularization scheme must enter into 
this computation but our conclusions depend explicitly 
from a cut-off A. Firstly, we point out the the sums on n,- 
remove some singular behavior being zero in such cases 
and the explicit physical cut-off A, in the infrared limit, 
is strictly linked to the string tension that is obtained 
from experiments. So, the argument for confinement 
appears consistent at this stage. 

5. Conclusions 

We provided a strong coupling expansion both for clas- 
sical and quantum field theory of a massless quartic 
scalar field and pure Yang-Mills theory. A set of classi- 
cal solutions is proved to exist for the Yang-Mills field, 
instantons, that support the view of a trivial infrared 
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